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Photoevaporation of protoplanetary discs II: evolutionary models 
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ABSTRACT 

We present a new model for protoplanetary disc evolution. This model combines viscous evo- 
lution with photoevaporation of the disc, in a manner similar to Clarke, Gendrin & Sotomavo^ 
However in a companion paper (Alexander. Clarke & Pringle 2006a) we have shown 
that at late times such models must consider the effect of stellar radiation directly incident on 
the inner disc edge, and here we model the observational implications of this process. We find 
that the entire disc is dispersed on a time-scale of order lO^yr after a disc lifetime of a few 
Myr, consistent with observations of T Tauri (TT) stars. We use a simple prescription to model 
the spectral energy distribution of the evolving disc, and demonstrate that the model is consis- 
tent with observational data across a wide range of wavelengths. We note also that the model 
predicts a short "inner hole" phase in the evolution of all TT discs, and make predictions for 
future observations at mid-infrared and milUmetre wavelengths. 

Key words: accretion, accretion discs - circumstellar matter - planetary systems: protoplan- 
etary discs - stars: pre-main-sequence 



1 INTRODUCTION 

For over twenty years the evolution, and eventual dispersal, of 
, discs around young stars has been an important area of study, 
' for theories of both star and planet formation. It is now well- 
established that at an age of ~ 10''yr most stars are surrounded by 
discs that are optically thick at optical and inf rared wavelengths 
iStrom et aDll989t iKenvon & HartmamJll995h . Observations at 
millimetre wavelengths show that these discs have masses that are 
typically a few percent of a solar mass iBeckwith et al. 1990), and 
so discs are widely believed to be potential sites for planet forma- 
tion. However at an age of ~ lO'yr most stars are not seen to have 
discs, suggesting t hat disc lifetimes are typically a few Myr (e.g. 
iHaisch et alj200lh . How stars lose their discs remains an unsolved 
question. 

The distribution of T Tauri stars (TTs) at infrared wavelengths 
provides some insight into this problem. These objects tend to fall 
into two distinct groups: those whose emission is consistent with 
a stellar photosphere plus an optically thick disc, and those which 
are compatible with purely photospheric emission. [Note that these 
two classes usually coincide with the spectroscopic classifications 
of classical and weak-lined T Tauri stars (CTTs and WTTs) respec- 
tively.] A number of authors have noted that very few transition ob- 
jects are observed between the CTT and WTT loci ( Skrutskie e^lj 
[l990: Kenvon & Hartmann 1995: Armitage, Clarke & Toul. 19991 
iH^taiam^^nilOOSh . Observations at longer wavelengths show 
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a similar behaviour, both in the mid-infrared jPersi et all 120001 
'Bontemos et al. 2001) and at millimetre wavelengths iDuvert et alJ 
2OO0: Andrews & Williams 2005). These observations suggest that 
discs are di spersed very rapidly, with the dispersal time estimated to 
be ~ lOVr isimon & Pratoll995l:lm)lk & Waltenl996) . Moreover 
the simultaneous decline is disc emission across such a wide range 
in wavelength suggests that the dispersal is essentially simultane- 
ous across the entire radial ext ent of the disc (see also discussion in 
iTakeuchi. Clarke & Litj2005h . 

This "two-time-scale" behaviour is inconsistent with conven- 
tion al models of dis c evolution, such as viscous evolution mod- 
els Jltotmam^^lJQ998) or models of magnetospheric clear- 
ing I Armitage et alJll999iV Such models predict power-law de- 
clines in disc properties, and therefore predict dispersal times 
whic h are always of the same order a s the disc lifetime. However 
I Clarke. Gendrin & SotomavoJ EooH hereafter CGSOl) showed 
that models which combine photoevaporation of the disc with 
viscous evolution can reproduce this two-time-scale behaviour. 
In this model, known as the "UV-switch" model (CGSOl), ion- 
izing radiation from the central star produces a photoevapo- 
rative wind at large radii. Detailed mo dels of photo evapora- 
tive winds were constructed by Hollenbac h et alj i 19941 see also 
iHollenbach. Yorke & Johnstoneir200Q) for the cases of both weak 
and strong stellar winds. In the case of TTs we consider only the 
weak stellar wind case. In this case ionizing radiation from the 
star creates an ionized layer on the surface of the disc, with con- 
ditions akin to an H ii region. Beyond some critical radius, known 
as the gravitational radius, the local thermal energy of the ionized is 
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greater than its gravitational energy and the gas escapes as a wind. 
The gravitational radius is therefore given by 



GM, I M. \ 



(1) 



where Cs is the sound speed of the ionized gas, typically lOkms 



iHollenbach et alj il994) show that the wind rate is determined by 
the the density at the ionization front, and find an integrated mass- 
loss rate of 
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where fl) is the ionizing flux produced by the star. More re- 
cent studies have extended the study of t he details in a number 
of way s ([Richling & Yorke 1997; Johnst one. HoUenbach & Ballvl 
ll998tlFbnt et aL.2004i') . and recent hydrodynamic modelling has 
resulted in slight modification of the quantitative results. When hy- 
drodynamic effects are considered the "efi'ective R^" is reduced by a 
factor of 5 (Liffman 2003; Font et a l. 2004J, and the mass-loss rate 
is reduced by a factor of around 3 |^onre^ij200^;_^foweve£_the 
qualita tive behaviour is unchanged from that of IHollenbach et alJ 

The so-called "UV-switch" model of CGSOl couples a photo- 
evaporative wind to a simple disc evolution model. At early times 
in the evolution the accretion rate through the disc is much larger 
than the wind rate, and the wind has a negligible effect. However 
at late times photoevaporation becomes important, depriving the 
disc of resupply inside Rg. At this point the inner disc drains on 
its own, short, viscous time-scale, giving a dispersal time much 
shorter than the disc lifetime . A number of similar studies have 
now been conducted ^Matsuvama, Johnstone & Hartmann 
| 2003l [Armitage, Clarke & Palla 2003; Ruden 2004; 
n^euchi. Clarke & LinI l2005t) . and this class of models show 
a number of attractive properties. 

However CGSOl highlighted two key problems with the 
model. Firstly, the model requires that TTs produce a rather large 
ionizing flux, of order 10*' ionizing photons per second. They also 
found that the outer disc, beyond was dispersed much too slowly 
to satisfy millimetre observations of WTTs, a finding re-affirmed by 
recent sub-millimetre observations I Andrews & Williams 2005). 
We have previously shown that it is reasonable to treat TT 
chromospheres as having a constant ionizing flux in the range 
~ 10'"-10'"*s-' jAlexander. Clarke & Pringi3i2005h . and we now 
seek to address the "outer disc problem". In a companion paper 
I Alexander, Clarke & Pringle 2006a, hereafter Paper I), we high- 
lighted an important flaw in the UV-switch mod el. The UV-switch 
model relies on the wind parametrization of IHollenbach et alj 
il994h . which assumes that the disc is extremely optically thick 
to Lyman continuum photons at all radii. [HoUenbach et alj il994l) 
find that the diffuse (recombination) field dominates the photoevap- 
oration at all radii of interest, as the direct field suffers extremely 
strong attenuation by the disc atmosphere. However we note that 
at late times in the UV-switch model the inner disc is drained, and 
is therefore optically thin to ionizing radiation. Consequently the 
direct field is important after the inner disc has drained. In Paper I 
we constructed detailed hydrodynamic models of the wind driven 
by the direct field, and derived a functional form for the mass-loss 
rate: 
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Here C and D are order-of-unity scaling constants, /i is the mean 
molecular weight of the gas, a is a power-law index, H/R is the 
ratio of the disc scale-height to radius, and Ri„ and Rout are the in- 
ner and outer disc radii respectively. Our numerical analysis fixed 
the values of the scaling constants to be a = 2.42 ± 0.09 and 
(CD) = 0.235±0.02 (for H/R = 0.05). As noted in Paper I, both the 
geometry of the radiative transfer problem and the form of the wind 
are qualitatively similar to the strong wind case of HoUenbach et 
{ 1994). However the effect of the direct radiation field is to increase 
the efficiency of the wind. Consequently the mass-loss rate due to 
direct photoevaporation is around an order of magnitude larger than 
that from the diffuse field, and is significant at late stages of the 
evolution. 

The diffusion equation for the evolution o f disc surface 
density E(i?, /) is j Lvnden-Bell & Pringld 1 1 9741: IPringld Il98ll: 
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where v is the kinematic viscosity and the term t,„i„^(R,t) repre- 
sents the mass-loss due to photoevaporatio n. CGSOl solved thi s 
equation using the "weak-wind" profile of HoUen bach et aljjl994b . 
At some point in the evolution the mass-loss rate from the wind 
falls to a level comparable to the accretion rate through the disc, 
and at this point the disc is rapidly drained inside R^. However, as 
mentioned above, CGSOl neglect the influence of the direct radi- 
ation field after this inner draining occurs. Consequently they find 
that the time-scale for dispersal of the outer disc is limited by the 
time material takes to diffuse inward to Rg (as the most of the mass- 
loss occurs close to Rg). Further, the R^^'^ dependence of the wind 
profile at large radii means that the mass-loss rate due to photoe- 
vaporation decreases significantly with time as the inner edge of 
the disc moves outward. Consequently the dispersal of the outer 
disc occurs on the viscous time-scale of the outer disc, and thus the 
outer disc is dispersed in a time comparable to the disc lifetime, 
much too slowly to satisfy observational constraints. We suggest 
that photoevaporation by the direct radiation field results in a dis- 
persal time significantly shorter than that predicted by CGSOl, and 
now seek to model the effects of this process on the evolution of the 
outer disc. 

In this paper we seek to incorporate the result of Paper I into 
models of disc evolution. We do this by solving the equation for 
the evolution of the surface density of a geometrically thin disc, 
including the photoevaporative wind as a sink term. In Section |5| 
we present a simple time-scale analysis, which demonstrates the 
significance of the direct radiation field. In Section^we construct 
a numerical model of disc evolution in the presence of a photoe- 
vaporative wind. We the use a simple prescription to model the ob- 
served spectral energy distribution (SED) of the evolving disc, and 
construct a set of models which cover a broad range in parameter 
space (Section|4j. In Section|5|we present our results, comparing 
the predicted SEDs to recent observational data. In Section |5| we 
discuss the implications and limitations of our results, and in Sec- 
tion0we summarize our conclusions. 



2 EVOLUTIONARY TIME-SCALES 

It is useful at this point to consider the evolutionary time-scales 
predicted by the different models. In the original model of CGSOl, 
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the inner disc drains on approximately the viscous time-scale at 
the draining radius, ty(Rg), which is significantly shorter than the 
disc lifetime to that point. However the outer disc drains on a much 
longer time-scale, comparable to the disc lifetime. After the inner 
disc has drained there is no accretion on to the star, so the only 
mass-loss is due to the wind. Consequently, if we neglect viscous 
evolution during clearing, the time-scale for the wind to clear the 
disc out to a radius R > Rg, t^^R), is given by 



Mj(< R) 



(5) 



where Md(< R) is the total disc mass at radii < R and Mwind(^) 
is the total mass-loss rate from the wind, integrated from radius R 
outwards. Neglecting numerical factors of order unity, the viscous 
time-scale for the evolution of an accretion disc is given by 



Mi(< R) 
Mi 



(6) 



where Mj is the disc accretion rate. We can therefore combine 
Equations|5|&|6|and express the clearing time as 



t,{R) = t,{R) 



Mwi„d(i?) 



(7) 



Clearing occurs when the disc accretion rate Mj falls to a level com- 
parable to the diffuse wind rate (which is constant, see Equation|2j, 
and therefore the clearing time-scale depends only on M„i„(j(i?). If 
^wind(^) is a decreasing function of R then accretion will dominate 
over the wind as the disc evolves, whereas if M„[ni{R) increases 
with R then the wind will dominate the evolution. 

In the model of CGSOl, where the diffuse field mass-loss 
profile is adopted throughout, the wind profile takes the form 
2wind(^) (R/R„y^^^. Consequently the total wind mass-loss rate, 
integrated from an inner edge radius R outwards, is My,ir,6(R) «^ 
if"''^, and we can normalise by noting that t^Rg) = ty(Rg). There- 
fore the clearing time-scale is given by 

>l/2 



t,(R) = ty(R)\^ 



(8) 



Thus we see that in the UV-switch model (CGSOl) the time-scale to 
clear the disc out to a radius R is longer than the viscous time-scale 
for all R > Rg. Consequently in this model viscosity dominates: 
mass is lost only after it has had time to diffuse inwards towards 
the inner disc edge. Thus while the inner disc satisfies the two- 
time-scale behaviour demanded by observations (see Section 0, 
the outer disc is dispersed much too slowly to satisfy the data. 

However if we consider direct photoevaporation of the outer 
disc we see a very different behaviour. In this case, as seen in Equa- 
tion|3|(see also Paper I), the integrated mass-loss from a disc with 
an inner edge at radius R is M„i„i oc R'-/^. Consequently in this case 
the clearing time-scale is given by 



tAR) = ty(R) 



-1/2 
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Thus we see that in the case of direct photoevaporation viscosity 
becomes progressively less significant as the inner edge evolves 
outwards, and that the evolution is instead dominated by the wind. 
Consequently we predict a much faster dispersal of the outer disc 
than that originally predicted by CGSOl. For TT parameters the 
draining radius is of order AU, and TT discs are observed to be up 
to several hundred AU in size. Thus we predict a clearing time for 
the entire disc that is ~1~10% of the viscous scaling time. Such 



behaviour satisfies the two-time-scale constraint across the entire 
radial extent of the disc. 



3 DISC MODEL 

We solve the diffusion equation for the disc surface density 
(Equation using a standard first-order explicit scheme (e.g. 
IPringle. Ver bunt & Wade 1986), w ith a g rid of points equispaced 
in i?''^. Following iHartmann et all il998h and CGSOl, we adopt a 
kinematic viscosity v which scales linearly with radius, so 



v(R) = vo 



R 
R'o 



(10) 



for some scaling value vo(/?o) at a scale radius Rq. Such a viscos- 
ity la w is consistent with an ff-prescription I Shakura & Sunvae^ 
1 19731) if the disc temperature at the midplane scales as 
and arguments in fav our of this viscosity law are discussed in 
IHartmann et alj il998h . We adopt an initial surface density pro- 
file consistent with the similarity solution of the diffusion equa- 
tion (Lvnden-Bell & Prinale 1974; Hartmann et al. 1998). Again 
following CGSOl, we define this profile to have the form 



I.(R): 



Mm 
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for an initial disc mass Md(0). In this form 1/e of the disc mass is 
initially at R > Rq, with an exponential decline in surface density 
at radii beyond Rq. However in practice the results are not sensitive 
to the form of the initial profile, and this form is chosen primarily 
as a simple means of parametrizing the initial disc mass. With this 
viscosity law the viscous scaling time, which governs the evolution 
of the disc, is given by 



3vo 



(12) 



and consequently the initial accretion rate at the origin is given by 



Mi(R = 0,1 = 0). 



3Md(0)vo 
2Rl 



(13) 



This fixes the scaling constant vq, and thus the disc model is entirely 
specified by the three parameters Md(0), Md(0, 0) and ^o- Initially 
we adopt Md(0) = O.OSMq, Md(0, 0) = 5.0 x lO^^Mgyr-' and Rq = 
lOAU. We use 1000 grid points, equispaced in i?''^, which span 
the radial range [0.0025 AU, 2500AU]. We also adopt adopt zero- 
torque boundary conditions throughout (i.e. we set E = at the grid 
boundaries) but note that the spatial domain is always large enough 
that the outer boundary condition has no effect on the results. 

At early tim es in the evolution we use the wind profile of 
iFont et al] i2004l kindly provided in numerical form by Ian Mc- 
Carthy), which incorporates more detailed hydrodynamics than the 
original work of Holl enbach et alJ iil994l) . However once the inner 
disc has been drained we alter the mass-loss profile to reflect the 
influence of direct photoevaporation on the outer disc. Thus it is 
necessary to define a numerical criterion for "draining", in order to 
indicate when to change to the direct profile. Additionally, as the 
form of the direct profile is normalised at the inner disc edge it is 
necessary to define the inner edge numerically also. 

By inspection of the mass-loss profiles we see that in both the 
diffuse and direct cases the mass-loss rate scales as 



M oc O 
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Evolution of surface density: M_^— IMq, 0—10 ' 




Figure 1. Evolution of surface density in model which incorporates direct 
photoevaporation. Snapshots of the surface density are plotted at ? = 0, 
2.0, 4.0, 5.9, 6.0, 6.01, 6.02, 6.02. . .6.18Myr. At t = 6.20Myr the surface 
density is zero across the entire grid. After the inner disc is drained direct 
photoevaporation disperses the outer disc very rapidly. 



Numerical analysis of the scaling constants shows that for equal 
ionizing fluxes the direct mass-loss rate exceeds that due to the dif- 
fuse field by a factor of 8.8 (for M, = IMq). Thus if the direct 
flux reaching the disc is greater than approximately 0.01 of its total 
value the direct wind exceeds the diffuse wind. Consequently our 
criterion for the "transition" between the two wind profiles is that 
the optical depth to ionizing photons, t, along the disc midplane to 
the draining radius (0.2/?g) satisfies 



exp(-T) = 0.01 



: 4.61 , 



(15) 



We evaluate the optical depth as r = Ncru.tev, where A' is the col- 
umn density along the disc midplane. Additionally, once the inner 
disc has drained we define the inner edge Ri„ to be the radius at 
which the optical depth reaches this critical value. After draining 
we adopt the mass-loss profile derived in Paper I (for H/R = 0.05): 



2„ind(^) = 2CDjjmHC,nin I — 

\«in 



where 



(16) 



(17) 



Here Ug is the Case B recombination coefficient for atomic hy- 
drogenat 10*K, which has a value of = 2.6 x 10^'^^cm^^s"' 
We adopt the best-fitting scaling constants from Pa- 
per I of (CD) = 0.235 and a = 2.42, and set ^ = 1.35 (following 
iHoUenbachet all 1994 CGSOl). 



3.1 Results 



same disc model was run with H/R = 0. 1 (and therefore best-fitting 
parameters of (CD) = 0.60 and a = 4.50, see Paper I). In this case 
the inner disc draining is identical to the case of H/R = 0.05, as the 
early evolution depends only on the diffuse wind and is independent 
of H/R. The outer disc is cleared somewhat more slowly by the 
direct wind than for H/R = 0.05, with the outer grid radius not 
reached until t = 6.27Myr. However this represents a change of 
less than a factor of two in the time required to clear the outer disc, 
so we are satisfied that the choice of H/R is not a significant factor 
in evolution of the outer disc. 

A second version of this model was run with v oc in 
order to investigate the effect of varying the viscosity law. The nu- 
merical scaling const ants in the similarity solut ion vary with th e 
viscosity law (Lvnden-Bell & Prinele..l974:.Hartmann et al^l998^ ■ 
but while the details differ the qualitative behaviour of the model 
is unchanged. The steeper surface density profile that results from 
this new viscosity law means that the inner disc draining is some- 
what slower than in the reference model, but this is countered by a 
more rapid clearing of the outer disc (as a smaller fraction of the 
disc mass now resides at large radii). The disc is still dispersed on a 
time-scale some 1-2 orders of magnitude shorter than the disc life- 
time', and the only significant consequence of changing the viscos- 
ity law is the expected modification of the surface density profile. 

Thus, as predicted in Section|2| when direct photoevaporation 
is taken into account the entire disc is dispersed on a time-scale 
approximately 2 orders of magnitude shorter than the disc lifetime. 
Thus it seems that this model has solved the "outer disc problem" 
that affected the original UV-switch model of CGSOl. In order to 
compare to observed data, however, it is necessary to model the 
observable properties of the disc. 



4 OBSERVABLE CONSEQUENCES: SPECTRAL 
ENERGY DISTRIBUTIONS 

In order to compare the results of the model to observed data it 
is necessary to use the model to generate spectral energy distri- 
butions (henceforth SEDs). Following Hartmann et al. 1 1998) and 
CGSOl, we assume that the disc is vertically isothermal and emits 
as a blackbody. Therefore the flux emitted by the disc at frequency 



COS i r 



InRBy (T(R)) 1 - exp 



'(- — )1 

\ cos I Jl 



dR, 



(18) 



where d is the Earth-star distance (taken to be 140pc, the distance 
of the Taurus-Auriga cloud), By(T) is the Planck function and i is 
the inclination angle of the disc (i = is face-on). T(R) is the radial 
temperature profile of the disc, and the optical depth Ty is evaluated 



T,, = K,I.(R) . 



(19) 



We adopt a st andard power-law for the dust opacity Ky 
jBeckwithetalll99m : 



FigQ shows the evolution of the surface density for the reference 
model. This model has M. = IMq, Md(0) = O.OSMq, Md(0,0) = 
5.0 X 10^'Moyr"' and /?o = lOAU. The results shown are for O = 
10''^s~'. The inner disc begins to drain at f = 6.01Myr, and the disc 
is cleared to the outer grid radius (2500AU) at t = 6.20Myr. At 
the point where the inner disc begins to drain the total disc mass 
remaining is approximately 0.002Mo, or approximately 2Mj„p. 
In order to investigate the dependence on the value of H/R the 



Ky = 0.1 



lO'^Hz 



2 -1 

cm^ g ' . 



(20) 



' Here the term "lifetime" refers to the age of the disc, t, when the disc is 
dispersed. 

^ Note that we have now used the symbol v to denote both frequency and 
kinematic viscosity, v denotes frequency only in Eg nations 1 1 8II20I and it 
should be obvious from context to which quantity the symbol refers. 
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We add the stellar contribution to the SED as the total flux emit- 
ted by a blackbody of temperature T, and radius R,. The values of 
these two parameters are taken from the pre-main sequence mod- 
els of Fibut et al. (1999, kindly provided in electronic form by 
Chris Tout). TT stars are contracting due to gravitational collapse, 
so for a given stellar mass both T, and R, vary somewhat with 
age. We adopt the median age of the T aurus-Auriga cloud, 2Myr 
jPalla & Stabler 2000: Hartmann'2001^■ throughout. 

Obviously the disc temperature adopted is crucial to the re- 
sultant SED. Further, we wish to investigate the effect of varying 
the stellar mass in our models, so it is necessary to define T(R) in 
such a way that it only dep ends on the stellar mass M,. We adopt a 
"flared reprocessing disc" jKenvon & Hartmannlll987l) power-law 
profile: 

T{R) = Tol — \ , i?>/?D (21) 

where the normalisation condition is set by the dust destruction ra- 
dius Rq. At radii smaller than this the temperature is set to zero, as 
no dust can survive here and we assume that the opacity due to the 
gas alone is negligible. We adopt a dust destruction temperature of 
Tj) = 1500K. Additionally, we adopt a minimum disc temperature 
of lOK to accoun t for external heating o f the disc (e.g. cosmic rays 
or diffuse UV, see'Hartmann et a\'l99^. By assuming that the stel- 
lar irradiation balances the emission from the vertical disc edge at 
Ro we find that 

i?D=A/;.|^J , (22) 

where A is a constant of proportionality that reflects how efficiently 
the disc radiates. 100% efficiency (the "small grains approxima- 
tion") gives A = V2, but by comparison to the observed "median 
SED" for Taurus- Auriga Jd'Alessio et alJlil999. see Fig|2| below) 
we find that A = 1.75 provides a better fit to the data. In adopting 
a value of A > V2 we essentially assume that the disc edge at Ro 
does not radiate as a blackbody (i.e. it radiates with less than 100% 
efficiency), an assumption verified by detailed radiative transfer 
models (e.e. ld'Alessio et all2005lJ) . We note at this point that the 
disc temperature at the "draining radius" of ^ 0.2Rg does not vary 
significantly with stellar mass. Obviously scales linearly with 
M,, but the effect of this on the disc temperature is offset by the 
variation of T, with stellar mass. We find that T(0.2Rg) = 540K for 
M = 0.2Mo, but is only slightly smaller (410K) for M = l.OMg. 
Consequently, while the draining radius scales linearly with stellar 
mass, the colour change resulting from inner disc draining is not 
especially sensitive to stellar mass. 

In addition to evaluating the SEDs, we use the predicted SEDs 
to generate magnitudes in the various photometric bands of the 2- 
Micron All Sky Siir\>ey (henceforth 2MASS) and the Spitzer Space 
Telescope. We use the filter transmission functions given on the 
2MASS and Spitzer Science Center websites^^, and use the given 
zero-point fluxes to convert these fluxes to magnitudes. 

In order to test this model we compare the predicted SED 
with the "median SED" of CTTs in Taurus-Auriga ( d' Alessio et al. 
11999 ). We adopt stellar parameters of T, = 3200K and R, = 1 .5Ro, 
consistent with M, = 0.3~0.4Mo (depending on age), and use the 
surface density profile of the reference disc at / = 0. As seen in 



^ |http : //www ■ ipac ■ caltech . e(iu/2mass/| and 
Pr»^77ss^^p^ze^^a^ecl^3u7 ^respectively. 
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Figure 2. SED produced by our model at f = for T, = 320OK and 
Rt = I.5R0, values consistent with a pre-main-sequence star of mass 0.3- 
O.4M0. The SED produced by the model is shown as a solid line, with the 
individual contributions from the stellar black-body and the disc shown as 
dashed lines. The points an d error bars are the "m edian SED" of CTTs in 
Taurus-Auriga, taken from ld'Alessio et all iiggsft . The two points at the 
shortest wavelengths are for U and B band data, and show the UV excess 
typical of CTTs. 

Fig|2| the model reproduces the observed data well out to a wave- 
length of ^ lOOyum, but rather under-predicts the flux at millime- 
tre wavelengths. This apparent error is not of great concern, how- 
ever, as in the initial disc configuration mass is mostly confined to 
small radii, with the result that the disc is optically thick at mil- 
limetre wavelengths. At later times the disc expands and becomes 
optically thin, boosting the millimetre flux. Additionally, once the 
disc is optically thin the emitted flux is very sensitive to the to- 
tal disc mass, and our reference model has rather a low initial 
disc mass (O.OSMq). More massive discs at later times (once the 
disc has spread beyond the scale radius Rq) match observed mil- 
limetre fluxes much better. In the Spitzer IRAC bands, our model 
gives colours of [3.6] - [4.5] = 0.46, [4.5] - [5.8] = 0.57, and 
[5.8] - [8.0] = 0.91. These compare favourably with the values 
for the median SED ([ 3 6] - [4.5] = 40, [4 .5] - [5.8] = 0.52, 
and [5.8] - [8.0] = 0.83. lHartmann et all2o'oa) . suggesting that the 
model predicts IRAC colours to an accuracy of around ±0.1 mag. 

4.1 Model set 

In order to study the behaviour of the SED a series of disc models 
were run with different stellar masses. As seen above, the outer disc 
evolution is not especially sensitive to the value of HIR, so when 
considering direct photoevaporation we adopt the best-fitting wind 
profile for HIR = 0.05. There are therefore five free parameters in 
the disc model: stellar mass M,, initial disc mass M^iO), initial ac- 
cretion rate M^iO, 0), scaling radius Rq and ionizing flux <t>. There- 
fore it is necessary to evaluate these parameters self-consistently in 
order to study the effect of a single parameter. In order to achieve 
this we assume that both the initial disc mass and the disc scal- 
ing radius scale linearly with M,, and we adopt the normalisation 
conditions 

/.o=10AU^, (23) 
and 
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Md(0) = 0.15M. , 



(24) 



The initial accretion rate is dependent on the disc viscosity. As 
noted above, we adopt a viscosity law of the form v oc R. We now 
fix the normalisation by assuming that that the viscosity parameter 
a Ishakura & Sunvaev 1973) and disc aspect ratio HjR are inde- 
pendent of M,. The orbital and viscous time-scales at are related 
by 

The relationships between the accretion rate, viscosity and viscous 
time-scale were given in Equations 1121 & [T3l We combine these 
expressions to fix the initial accretion rate as 



Md(0,0) 



Md(0) 

2ty 



(26) 



For typical parameters EQuation l25l gives ty =^ lOOO/oit, where is 
simply the Keplerian orbital time at Rg. (This gives f,, = 3.2 x 10''yr 
at Ro = lOAU.) Thus the disc parameters are specified in a man- 
ner which depends only on the stellar mass M,, due to the manner 
in which Md(0) and Rq are specified. This fiducial disc model was 
evaluated for stellar masses of M, = 0.2, 0.5, 1.0 & 2.0Mo, and 
for each model the SED was evaluated as a function of time at an 
inclination angle of ; = 60° (i.e. cos i = 0.5, the mean inclination 
of a random sample). We set O = lO^'^s"' in all of these mod- 
els. In order to explore parameter space further four additional disc 
models were run for M, = I.OMq, with Rq = 5AU, O = 10*'s"', 
Md(0) = 0.15M. and = 5000/orb- Lastly the SED for the fiducial 
IMq model was also evaluated for an inclination angles of ! = 
(i.e. a face-on disc) and / = 80° (i.e. a nearly edge-on disc). Thus 
one group of models explores the effect of stellar mass on the ob- 
served disc emission, while the second group explores the conse- 
quences of varying the disc parameters for a fixed stellar mass. 

Magnitudes were evaluated in the J, H and 2MASS bands, 
and the four IRAC bands (which have central wavelengths of 3.6, 
4.5, 5.8 & 8.0yum respectively). In addition magnitudes were cal- 
culated in the Spitzer 24/jm MIPS band''. The flux at 850yum was 
also measured, for comparison to sub-millimetre (SCUBA) ob- 
servations. Each of these four wavebands (observed by 2MASS, 
IRAC, MIPS and SCUBA) probes a different region of the disc, 
with each providing different constraints. The 2MASS bands pri- 
marily observe the stellar flux, and therefore provide a valuable 
normalisation condition. The IRAC bands probe the inner disc, at 
the draining radius and smaller radii, while the 24/jm MIPS band 
probes the emission at somewhat larger radii, beyond the initial 
draining radius. The emission in both the IRAC and MIPS bands is 
(mostly) optically thick (the longer wavelength bands show weak 
optical depth effects), and so is sensitive only to the disc temper- 
ature and to whether or not the disc has drained. Therefore small 
changes in the disc model do not have a significant effect here. 
Lastly the flux at 850/jm, as measured by SCUBA, measures op- 
tically thin emission from the entire disc, and is therefore rather 
sensitive to a number of the parameters in the disc model. 

We now seek to compare the results from our disc models 
to observed data. In order to do this we have created a compos- 
ite dataset, using observations taken from the literature. We use 



Note that the zero-point flux in the 24/.im MIPS band, 7.3Jy, is rather 
uncertain. Consequently the absolute value of the magnitudes in this band 
are subject to systematic errors. However this merely shifts the zero-point 
of the magnitude system, and any trends are unaffected. 
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AU 
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1.0 


1.0 


0.15 


10.0 


1000 


8.37 


8.47 


0.2 


1.0 


0.15 


2.0 


1000 


3.52 


3.57 


0.5 


1.0 


0.15 


5.0 


1000 


6.13 


6.21 


2.0 


1.0 


0.15 


20.0 


1000 


10.66 


10.79 


1.0 


10.0 


0.15 


10.0 


1000 


4.87 


4.97 


1.0 


1.0 


0.3 


10.0 


1000 


10.56 


10.67 


1.0 


1.0 


0.15 


5.0 


1000 


5.91 


5.98 


1.0 


1.0 


0.15 


10.0 


5000 


18.47 


18.96 



Table 1. Table showing the parameters used in the various disc evolution 
models, and also the disc lifetimes found for each model. The column la- 
belled fi indicates the time at which the inner disc was drained (i.e. when the 
model switches from the diffuse to direct wind profile), while the column 
labelled f2 indicates the time at which the entire disc was dispersed. 



data from recent sur veys of the Taurus-Au riga cloud, from the 
IRAC observations of Hartmann et al. (2005) and the SCUBA ob- 
servations of Andrews & Williams 1 2005). I Hartmann et al. 2 00^ 
also list 2MASS magnitudes for all sources.) We include only 
the sources which are unambiguously included in both samples, 
rejecting any binaries which are resolved by IRAC but not by 
SCUBA. This leaves a total of 42 objects: 29 CTT s, 12 WTTs 
and t he possible transition o bject CoKu Tau/4 (see iForrest et alJ 
2004t Id' Alessio et aljEoOsS) . All but 3 of the CTTs are detected 
by SCUBA, with upper limits only found for DP Tau, CIDA 1 1 and 
CIDA 12. (These latter two objects are close to the brown dwarf 
limit and may not be "true" TTs.) By contrast only 1 of the WTTs 
(LkHQr332 Gl) is detected at 850;/m, confirming that, in general, 
disc dispersal occurs simu ltaneously over the entire ra dial extent of 
the disc (see discussion in lAndrews & Williamsl200^ . 



5 RESULTS 

The results of our models are shown in Table Q and Figs|3lj5] 
Table Q shows the disc lifetimes predicted by the model, which 
are entirely consistent with disc lifetimes of HOMyr and dis- 
persal times of order 10^ yr (as derived from observations, e.g. 
iKenvon & HartmannllT99a iHaisch et alj|200lh . Fig|3l shows evo- 
lutionary tracks on a - [3.6] / Ks - [8.0] two-colour diagram 
(analogous to previously published K — L I K - N plots, e.g. 
iKenvon & H artmann 199^. lArmitage et alj l999). The data points 
show a clear gap between the loci of CTTs and WTTs which the 
tracks reproduce well, showing a rapid transition across the gap. 
The disc emission is optically thick, and so we see Class II colours 
before the disc is cleared followed by Class III colours afterwards. 
We also see that stellar mass and disc inclination angle are the dom- 
inant effects at these wavelengths. For a given stellar mass and in- 
clination angle different disc models show very similar tracks, as 
the infrared emission is generally optically thick. Consequently the 
tracks depend primarily on the temperature profile, and are insensi- 
tive to the parameters of the disc model. Minor optical depth effects 
are seen at 8yum in the models with O = 10*'' s"' and ty = 5000fo,b! 
but these are not significant. The models struggle to reproduce the 
extreme points in the CTT distribution, at both the red and blue 
ends, and indeed the A^s-[8.0] colours of our model during the CTT 
phase are rather redder than the observed data. However there are 
several factors which can account for this. The SED is rather sensi- 
tive to the adopted stellar temperatures, and the disc emission can 
also be increased by so-called "accretion luminosity". This arises 
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Figure 3. 2MASS/IRAC K^, - [3.6] / - [8.0] plots, with data points 
from Hartmann et al. i20()5). Solid circles represent CTTs and open circles 
WTTs, with the possible transition object CoKu Tau/4 represented by a star. 
The upper panel shows evolutionary tracks for different stellar masses with 
inchnation angle = 60°: M, = 0.2 (dotted hne), 0.5 (dashed), 1.0 (solid) 
and 2.OM0 (dot-dashed). The hght grey track is for M, = LOMq with = 0, 
and the dark grey track i = 80° . The lower panel shows the effect of varying 
the disc parameters with M, = I.OMq. The solid black track is the fiducial 
model (as in the upper panel). The M(j(0) = 0.3Mo track is shown in grey 
(obscured by the black track), <D = lO^'^s"' as a dashed line, Ro = 5AU as a 
dotted line, and ty = 5000forb as a dot-dashed line. In both plots points are 
added to the tracks every lO^yr to illustrate the evolution. 



due to viscous heating of the disc, and boosts the disc emission. 
It is omitted from our SED model but can be signi ficant in CTTs 
with high accretion rates fe.g. lArmitage e t al. 1999). Additionally, 
our model does not include emission from the inner disc edge (or 
"wall") at the dust destruction radi us, which may con tribute signif- 
icantly to the SED at /I ^ 3^m (e.g. lNatta etall200ll) . We note also 
that our models predict rather bluer Class III colours than those 
observed. This is due to the rather unrealistic use of a blackbody 
stellar spectrum, and we do not consider it to be a serious problem 
with the model. 

IVlore interesting is the evolution of the millimetre flux, as 
shown in Fig|31 The millimetre emission is mostly optically thin, 
and so is much more sensitive to the disc mass than the emission at 
shorter wavelengths. However as it is optically thin, it is relatively 
insensitive to the inclination angle. The behaviour of the tracks in 
the Fv(850yum) / - [8.0] plane, shown in Fig|4| is explained as 




K3-[8.0] 

Figu re 4. - [8.01 plotte d ag ainst 850/jm flux, with data points taken 
from lHartmann et alj l2005t) and lAndrews & Wilham'3 <2005t) . Filled sym- 
bols represent CTTs, open symbols WTTs, and the star CoKu Tau/4. Circles 
represent objects detected at 850/jm, while triangles denote 3<t upper hmits. 
As in Fig|3] the upper panel shows the effect of varying stellar mass in the 
models, while the lower panel shows the eft'ect of varying the disc param- 
eters. The line styles are the same as those used in Fig|3] with points again 
plotted every 10^ yr. 



follows. The initial evolution shows an increasing 850yum flux at 
fixed infrared colour, where the infrared colour is dependent on 
both stellar mass and disc inclination angle. The increase in 850/jm 
flux occurs because the disc viscously expands from its initial con- 
figuration, which is optically thick at 850yum. However soon ac- 
cretion takes over and the flux, now optically thin and simply pro- 
portional to the disc mass, declines at the total disc mass decreases. 
The tracks fall vertically (as the infrared emission is optically thick) 
until the inner disc is drained, at which point the Ks - [8.0] colour 
rapidly "jumps" to a stellar value at a fixed millimetre flux. The 
level of this flux is determined by the disc mass and temperature 
at this point, and does not change significantly as the inner disc 
drains. This is due to the v oc i? viscosity law, which forces most 
of the disc mass to reside at large radii. Draining occurs at a fixed 
value of the disc accretion rate, approximately equal to that of the 
diffuse wi nd. However the disc accretion rate can be expressed as 
M ~ vS (iPringlelll98ll) . Therefore for a fixed stellar mass we ex- 
pect that the level of the 850yum flux when the inner disc drains 
should depend only on the ionizing flux and the viscosity law. This 
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Figure 5. Predicted evolutionary tracks in tlie 2MASS/IRAC/MIPS ~ 
[4.5] I - [24] two-colour diagram. As in Fig|3]the upper panel shows 
the effect of varying stellar mass, while the lower panel shows the effect of 
varying the disc parameters. The line styles are again the same as in Fig|3l 
with points again plotted every 10^ yr. 



is verified in Fig|4| where the models with increased <5> and lower 
viscosity show significantly larger 850yum fluxes as the inner disc 
drains. (The total disc mass at this point in the fiducial model is 
Q.OQIMq. In the model with C) = 10*'s"' the disc mass at drain- 
ing is 0.006Mq.) At this point the evolution "stalls" for ~ lO^yr. 
This occurs because most of the disc mass resides at large radii, 
so the 850/im flux remains approximately constant while the inner 
disc is cleared. Once direct photoevaporation starts to clear the disc 
at large radii the 850pm flux falls rapidly to a very low level. The 
original UV-switch model (CGSOl), which omitted direct photo- 
evaporation, resulted in a millimetre flux of a few mJy which re- 
mained at very late times. This is clearly not a problem with our 
new model, which predicts millimetre fluxes of order 10^'mJy (the 
stellar contribution) once the disc has been cleared. However the 
model predicts that for around lO'yr objects should show stellar 
near- to mid-infrared colours but retain millimetre fluxes of 1- 
lOmJy. This may explain the two objects detected in this region 
by Andrews & Williams (2005): CoKu Tau/4 and LkHa332 Gl. 
These detections are near to the sensitivity limit of current obser- 
vations. Our model predicts that a factor of ~ 10 increase in sen- 
sitivity should results in the detection at millimetre wavelengths of 
an increased number of sources with Class III infrared SEDs, rep- 



resenting a few percent of the total population. This represents a 
valuable future test of this model. 

Another interesting result of the models is the predicted be- 
haviour in the mid-infrared at 20-50;um. This region of the spec- 
trum probes the region outside the draining radius, but in a wave- 
length range where the emission is (mostly) expected to be opti- 
cally thick. To date few observations have been made here, but this 
should be remedied in the near future by MIPS observations. Fig|5| 
shows the predicted tracks in the A", - [4.5] I K^- [24] two-colour 
diagram. Here we initially see a slow decline in the 24pm emis- 
sion due to weak optical depth effects. The "inside-out" manner 
of the disc clearing causes the emission in the near infrared to de- 
cline more rapidly than that at 24yum, with the models predicting 
that a significant population of sources, again at the few percent 
level, should show Class III colours in the IRAC bands but signifi- 
cant excesses at 24yum. Forthcoming observations from Spitzer will 
therefore provide another valuable test of the model. 

Lastly, we note that it is possible to constrain some disc pa- 
rameters from the model results. By demanding that disc li fetimes 
be in the 1-lOlVlyr range, as seen in observations (e.g. Hais ch et alJ 
'2001"), it is possible to constrain the disc viscosity, albeit rather 
weakly. For initial disc masses of 0.1 5M, the ratio of viscous to 
orbital times must be =^ 1000, and certainly less than 5000. For re- 
alistic H/R ratios this suggests a viscosity parameter in the range 
0.02 < c < 0.4. There is no real consensus as to what value of 
a can be produced by models of angular momentum tra nsport in 
discs, but fiducial values tend to be of order a ~ 0.01 tStone et alJ 
2000). Thus these derived values of a are near to the upper limit of 
those predicted by current models of angular momentum transport 
in discs, and may pose problems for these models. 

Similarly, by comparing the models to observed millimetre 
fluxes (as in Fig|3J it is possible to place indirect constraints on the 
ionizing flux O. The fact that few CTTs are not detected at 850yum 
at the lOmly level, and that similarly few WTTs are detected at 
the same level, suggests that the millimetre flux during the "transi- 
tion" phase must be in the 1-lOmJy range for most objects. (If this 
were not the case we would expect to see either many more WTTs 
detected or many more CTTs not detected at the lOmly level at 
850yum.) As seen in Fig|4|this level depends on both the disc and 
stellar parameters, and not solely on O. However the results suggest 
ionizing fluxes in the ra nge ~ 10'*'-10''^s~', a range consistent with 
our previous estimates iAlexander et alj2005l) . 



6 DISCUSSION 

There are obvious limitations to the model. Firstly, we note that 
the manner in which we treat the disc temperature and viscosity 
profiles are not self-consistent. The derived wind profiles (from Pa- 
per I) assume a constant H/R, but throughout this paper we adopt 
a disc temperature profile which implies disc flaring. Further, al- 
though the temperature profile adopted when evaluating the SED 
does result in a v oc ^ viscosity law, the normalisation of the two 
parametrizations is inconsistent. However the viscosity depends 
primarily on the midplane temperature, while the emitted SED de- 
pends on the surface temperature. T hese are most likely not the 
same fe.g. lChiang & Goldreichll997l) . and so such an inconsistency 
is not unreasonable. Therefore while there are minor inconsisten- 
cies in the disc model, they do have some physical motivation and 
they do not have a strong effect on the model results. 

In addition the SED model is rather simple, employing a 
power-law for the disc temperature structure and neglecting sev- 
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eral possibly important factors, such as viscous iieating or emission 
from tlie inner wall. Further, the temperature normalisation is rather 
sensitive to the stellar temperature adopted, as seen in EQuation l22l 
The flux emitted by the disc is very sensitive to disc temperature, 
so the output SEDs are in turn rather sensitive to both the stellar 
temperature and the normalisation condition. Our models struggle 
to reproduce the reddest and bluest CTTs in the data and, as seen in 
Fi^ our simple model predicts rather redder A's- [8.0] colours than 
those observed during the CTT phase. However this is easily reme- 
died by small alterations to the SED model, and we note that our 
simple treatment of the disc neglects several possibly important ef- 
fects (see Section|5j, such as accretion luminosity or emission from 
the inner disc wall. Additionally, our model uses a single temper- 
ature blackbody to model the stellar flux. Model atmospheres can 
difi^er markedly fro m blackbody spectra in the near-infrared (e.g. 
iBarafTe et al1ll99d) . and this may also affect the colours obtained 
from the SED model. The model also struggles to reproduce the 
handful of objects with the largest observed 850/jm fluxes. How- 
ever, as seen above, the emission at millimetre wavelengths is ex- 
tremely sensitive to the disc parameters. Additionally, these ob- 
jects tend to have rather unusual SED slopes in the sub-millimetre 
(Andrews & Williams 2005), and so we do not consider this to be 
a serious problem with the model. 

A further consideration is the manner in which the transition 
from the diffuse to direct regime is treated (see Sectionj^. We de- 
fine a critical value for the optical depth to ionizing photons along 
the disc midplane and use it to switch instantaneously between the 
two wind parametrizations. A more realistic treatment would grad- 
ually increase the strength of the direct field as the inner disc is 
drained. Additionally, the behaviour of the diffuse field as the inner 
disc drains should be considered in more detail, as the diffuse wind 
model assumes that the underlying disc is always optically thick to 
ionizing photons. The transition occurs very rapidly, and its treat- 
ment does not have a strong effect on the overall evolution of the 
disc. The transition is important, however, when we consider the 
detailed conseqences for the observed SED during clearing. This is 
because the inner disc becomes optically thin to infrared emission 
(perpendicular to the disc midplane) rather earlier than it becomes 
optically thin to Lyman continuum photons (along the midplane). 
Objects appear as "inner hole" sources during this phase, and the 
wavelength at which the disc emission "cuts in" is rather sensitive 
to the SED model. 

Recentlv lMcCabe et alj i2006h have observed a strong corre- 
lation in passive (non- accreting) discs, between stellar spectral type 
and the wavelength at which outer disc emission becomes visible. 
They define "passive discs" as sources which show excess emission 
at 10-12yum but lack any near-infrared excess, and in a survey of bi- 
nary secondary stars they find that that stars of later spectral type 
are more likely to possess passive discs than stars of earlier spec- 
tral type. The binary separations are large, so within the framework 
of a photoevaporation model they attribute this to the linear scal- 
ing of the draining radius with stellar mass. Thus low mass stars 
have smaller holes, which they argue should contribute disc emis- 
sion at shorter wavelengths than the corresponding holes in higher 
mass stars. This interpretation is predicated on the assumption that 
the smaller holes are hotter. However, in our models, smaller holes 
turn out to be only modestly hotter, since the stellar photospheric 
temperature (which sets the temperature of the reprocessing disc) 
is lower for less massive stars. Further, this interpretation remains 
valid only if the inner disc clearing takes significantly longer than 
the outer disc clearing. If the outer disc clearing dominates the dis- 
persal time then the distribution of observed hole sizes at any given 



time is independent of stellar mass. Therefore a more accurate treat- 
ment of the transition from the diffuse to direct wind is necessary in 
order to make such detailed predictions. We also stress that a better 
treatment of the disc SED, in particular the contribution from the 
inner wall, is needed before we can asse ss whether our models can 
reproduce the passive discs observed bv lMcCabeet"ai]<200dl . 

We note also that our model neglects the eff^ects of non- 
ionizing far-ultraviolet (FUV) radiation on the evolution of the 
disc. FUV radiation heats the disc to lower tem peratures than ion- 
izing radiation (^ lOOOK cf. ^ 10,000K, Jo hnstone et alJll998t 
iGorti & Hollenbacd 12004) and therefore FUV-driven winds are 
launched from rather larger radii than those driven by ionizing ra- 
diation, typically launching at radii from 20-lOOAU. Ada ms et alj 
(.2004.) have shown that external FUV radiation can produce a sig- 
nificant disc wind, which is is dominated by flow from the outer 
disc edge. However this wind relies on the presence of a strong 
interstellar FUV radiation field, which is only present in the imme- 
diate vicinity of massive O- or B-type stars. Self-consistent models 
of FUV heating from the central object have not yet been com- 
pleted, but it seems likely that FUV radiation can drive a wind 
from the outer disc edge in this case also. The interplay between 
mass-loss from the outer disc edge (from FUV heating), mass-loss 
from smaller radii (from Lyman continuum heating) and viscous 
evolution will provide an interesting area for future study. While 
the details of such a calculation are beyond the scope of this paper, 
we note that if the stellar FUV radiation field can drive a significant 
wind from the outer edge of the disc then this may significantly 
shorten the disc lifetime in any given disc model. Consequently we 
expect that the inclusion of such a wind will reduce the viscosity 
(a) required to satisfy observational constraints on disc lifetimes. 

Another important consideration is dust. Dust is responsible 
for all of the emission discussed above, and our SED modelling 
assumes a constant dust-to-gas ratio throughout. However it is not 
at all obvious that the dust in the disc remains coupled to the gas. 
There are several forces which act o n du st grains in discs, all of 
which may be important here (see e.g. lOusta fson 1994). In addition 
to gravity the dust grains feel a drag force from the gas, and grains 
can both grow and be ground down through coUisional processes. 
Grains can also feel a force due to radiation pressure. In most cases 
the optical depth through the disc is such that this force is negligi- 
ble, but near to the inner disc edge it may become significant. Ra- 
diation from the star is also responsible for the Poynting-Roberston 
(P-R) effect, which can cause grains to lose angular momentum 
and slowly spiral inwards. Again, this effect is not significant if the 
grains are shielded from the stellar radiation field, but may become 
significant once the gas disc begins to drain. Additionally, as noted 
in Paper I, any dust remaining in the "inner hole" may absorb some 
ionizing photons, reducing the efficiency of the wind and lengthen- 
ing the dispersal time-scale. 

The net effect of these forces on the dust in the disc is diffi- 
cult to predict without making detailed calculations. The evolution 
of dust in the original UV-switch model has been studied in some 
detail by Takeuchi et al. (2005). They find that gas drag causes 
millimetre-size dust grains to accrete on to the star more rapidly 
than the gas unless the grains are rather "fluffy" (i.e. have rather low 
densities). In this model the grains are removed from the outer disc 
before the inner disc is drained, and the decline in disc emission is 
primarily due to the rapid accretion of the grains rather than global 
evolution of the disc. Howev er the sub-millimetre observations by 
lAndrews & William j iiool do not detect a significant fraction of 
objects which have no dust in the outer disc but are still accret- 
ing gas from their inner disc, suggesting that the dust is removed 
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too quickly in the model of iTakeuchi et alj i2005h . Additionally, 
iTakeuchi et alJ (2005) find that the force exerted on the grains by 
the wind is always less than gravity, so little dust is carried away by 
the photoevaporative wind. This is not of great significance in their 
model, as the grains are removed from the outer disc by gas drag, 
but may be a problem when considering direct photoevaporation of 
the outer disc. In this case gas drag cannot migrate the grains be- 
yond the inner disc edge, and the wind cannot carry away any dust 
grains found in the outer disc. Consequently we would expect the 
dust grains in the outer disc to remain after the gas has been dis- 
persed. However the time-scales for grains to be removed by either 
coUisional processes or by P- R drag are very short, of order 10^^- 
lO'^yr at radii of a few AU (e.g. Gustafson 1994). Consequently any 
dust "left behind" in this manner is expected to be removed suffi- 
ciently quickly as to have a negligible efltect on the evolution of the 

SED. 

However the analysis of ITakeuchi et alj i2005l) does not con- 
sider the effects of radiation-powered forces on the grains. This is 
not a significant problem in their model as the dust is shielded from 
stellar irradiation, but may well become significant in our model 
once the inner disc beg ins to drain. lKlahr & LirJ<200lll2005h found 
that radiation pressure in discs with low opacity can result in a 
clumping instability in the dust grains. In our model this effect may 
well become important when the inner disc drains, as the optical 
depth is essentially zero inside the inner disc edge. Such clumping 
increases grain collision rates, and therefore may have important 
consequences for planet formation theories as well as on the ob- 
served SED. We also note that any dust exposed to the stellar ra- 
diation field in this manner will attenuate the direct radiation field 
responsible for the photoevaporative wind (see discussion in Pa- 
per I). Additional effects, such as grain growth and dust replen- 
ishment, may also be significant. In short, it is unlikely that the 
disc evolves with a constant gas-to-dust ratio, as we have assumed. 
However the manner in which the dust evolves is not at all obvious, 
and we make no attempt to model it here. We also note in pass- 
ing that a h andful of gas-poor discs have recently been observed by 
iNaiita & Williams 1 ,2005) , and suggest that photoevaporation may 
have played a significant role in their evolution. 

6.1 Inner holes 

Lastly, we note that one consequence of the model is that all ob- 
jects are predicted to pass through an "inner hole" phase, where 
the disc is drained close to the star but remains "normal" at larger 
radii. Several su ch object s have now been o bserved ( e.g. GM Aur, 
[ Rice et alfeOoS D M Tau.lCalvet et al J2005I : TW Hva. lCalvet etall 
l2002t CoKu'fauM. iForrest et all200 4'). and it has been suggested 
that these may represent a class of "transition objects", intermediate 
between the CTT and WTT states. It is not clear that these objects 
represent a homogeneous group (in fact this seems unlikely), but it 
is interesting to compare them to the predictions of our model. 

Despite the presence of holes in their (dust) discs, both GM 
Aur and DM Tau have l arge accretion rates (o n to the star), of or- 
der 10-*Moyr-' iCalvet e t al. 2005: IHueso & Guillot 2005). Ad- 
ditionally, in the case of GM Aur the obser ved inner hole i s in- 
compatible with a standard gas-to-dust ratio jRice et all20'03l) . As 
such GM Aur is obviously inconsistent with the photoevaporation 
model, and some other process, such as grain growth or the pres- 
ence of a planet, must be at work here. Both GM Aur and DM Tau 
have also been observ ed to have rel atively large outer disc masses, 
of order 0.05~0.IMo |c^vetetai]|2005 ; Hueso & Guillot 2005). 
In the photoevaporation model the disc mass at draining is approx- 



imately equal to the product of the (diffuse) wind mass-loss rate 
(~ lO^'^Moyr"') and the viscous evolution timescale (~ lO'yr), 
and, as seen in Section ITTI is typically O.OOlMg (plus or mi- 
nus around an order of magnitude). Therefore it is rather difficult 
(though not completely impossible) for photoevaporation to pro- 
duce inner holes with outer disc masses as large as those observed 
in these two sources, again suggesting that photoevaporation is not 
responsible for the holes observed in these two discs. However we 
note that in deriving these disc masses .Calvet et a l. ( 2005) have de- 
rived rather small viscosity parameters {a ~ 0.001), and that larger 
values of a would imply correspondingly smaller disc masses (see 
also Hueso & Guillot 2005). 

TW Hya has a "hole" in the disc inside 3~4AU, but still shows 
some dust and gas emission from inside this radius I Calvet et all 
200 2). It also has a me asured accretion rate of ^ 4 x lO^'^Mgyr"' 
iMuzeroUe et alJl2000l) . and it has been suggested that a planet 
may have formed (or be forming) in the disc. The disc mass is 
rather larger than that predicted by the photoevaporation model, but 
the low accretion rate is similar to that produced as the inner disc 
drains in the photoevaporation model. Therefore in principle it may 
be possible to model the SED with a photoevaporation model, al- 
though in order to reproduce the observed accretion rate the model 
requires that the object be observed during a rather short "window" 
in the evolution (while the inner disc is in the process of draining). 
Consequently it is not yet clear whether or not photoevaporation 
plays a role here. 

CoKu Tau/4, however, is almost entirely devoid of material 
inside lOAU and shows no evidence of a ccretion on to the star 
jporrest et alj |2004l: ^d' Alessio e t_aljr2005 j) . The outer disc mass 
of O.OOIMm jd'Alessio et al. 2001^ is similar to that predicted by 
our model (see Section ITTI and discussion above). We have seen 
above that the evolutionary tracks produced by the photoevapo- 
ration model can reproduce the SED of CoKu Tau/4 well, and at 
present this object seems to be entirely consistent with the predic- 
tions of the photoevaporation model. 

We note that much current work has invoked the presence of 
a formed or forming planet in order to explain such holes (e.g. 
ICalvet et alj200llRice et alj2003tlOuillen et alj2004 . Whilst this 
may well turn out to be the case, we emphasise that other mecha- 
nisms, such as photoevaporation, can also produce holes in discs, 
and that holes do not necessarily imply the presence of a planet. 
Indeed it seems likely that photoevaporation and planet formation 
occur simultaneously during the evolution of protoplanetary discs, 
and it is not yet clear whether photoevaporation acts to accelerate or 
retard planet formation. Spitzer is expected to discover many more 
objects similar to those discussed here, and only with more data 
will the nature of these objects become apparent. However the cur- 
rent data suggest that the inner hole sources discovered to date are 
far from a homogeneous class of objects. 



7 SUMMARY 

We have presented a new model for the evolution of protoplane- 
tary discs. Our model combines viscous evolution with photoevap- 
oration of the disc by stellar radiation, and includes the effects of 
the direct radiation field at late stages of the evolution (using the 
wind prescription derived in Paper I). We have constructed a nu- 
merical model for the evolving disc, and used a simple prescrip- 
tion to model the behaviour of the spectral energy distribution as 
the disc evolves. Our fiducial model predicts that the disc is com- 
pletely dispersed on a time-scale of order lO'yr after a lifetime of a 
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few Myr, consistent with observationally derived time-scales. Our 
results are consistent with observational data across a broad range 
in wavelengths, and allow us to place weak constraints on some 
model parameters: we derive ionizing fluxes in the range ~ 10"*'- 
lO^'^s"' and viscosity parameters (a) in the range ~ 0.02-0.4. We 
also make predictions as to what will be seen in future observations 
in the mid-infrared and at millimetre wavelengths. To date this is 
the only model of disc evolution which can reproduce the rapid 
disc dispersal seen in observations of T Tauri discs. We note also 
that the model suggests that all evolving discs pass through a short 
"iimer hole" phase. During this phase the outer disc mass is of order 
O.OOIM0, and we predict that inner hole sources should represent 
around 1-10% of the observed population of TTs. We compare the 
predictions of the model with the handful of inner hole sources ob- 
served to date and show that some, but not all, of these objects are 
consistent with the predictions of our model. 
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